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Abstract 

We present a formula for a classical r-matrix of an integrable system obtained by 
Hamiltonian reduction of some free field theories using pure gauge symmetries. The 
framework of the reduction is restricted only by the assumption that the respective 
gauge transformations are Lie group ones. Our formula is in terms of Dirac brackets, 
and some new observations on these brackets are made. We apply our method to derive 
a classical r-matrix for the elliptic Calogero-Moser system with spin starting from the 
Higgs bundle over an elliptic curve with marked points. In the paper we also derive 
a classical Feigin-Odesskii algebra by a Poisson reduction of some modification of the 
Higgs bundle over an elliptic curve. This allows us to include integrable lattice models 
in a Hitchin type construction. 

1 Introduction 



A classical r-matrix structure is an important tool for investigating integrable systems 
^miS]- It encodes the Hamiltonian structure of the Lax equation, provides the involution 
of integrals of motion, and gives a natural framework for quantizing integrable systems in 
a quantum group theoretic setting |3|. 

The aim of this paper is severalfold. First, we present a formula for the classical r- 
matrices of integrable systems, derived in the framework of Dirac's Hamiltonian reduction 
for systems involving only gauge transformations. In the process we shall derive new results 
describing the Dirac brackets. As an application of our general formula we shall calculate a 
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classical r-matrix for the elliptic Calogero-Moser system with spin. Following the works [SI 
"v^s start from the Higgs bundle over an elliptic curve with a marked point. This is a free 
theory with a trivial r-matrix. After the Dirac procedure we come to the desired r-matrix. 
Finally we derive the classical Feigin-Odesskii algebra [S] in the similar fashion. Here the 
corresponding unreduced space is similar to the cotangent bundle of the centrally extended 
loop group L{GL]\f). The Poisson structure on this space is a particular example of a general 
construction proposed in [0]. In this way we derive the corresponding quadratic Poisson 
algebra using the Hitchin approach ^U] . Thus we have managed to include integrable lattice 
models in the general Hitchin construction. 

Our work is motivated by the papers jlll I12j where the authors use a gauge invariant 
extension of Lax matrices to derive classical r-matrices for the Toda chain, and for the 
trigonometric and elliptic (spin zero) Calogero-Moser systems. In the first paper the 
Hamiltonian reduction is considered on the cotangent bundle to a finite dimensional Lie 
group, while in the second paper the construction is generalised to the case of the central 
extension of the loop group. In this context the paper also warrants mention. There the 
authors consider a special case of Poisson reduction on Poisson-Lie groupoids in order to 
obtain new examples of the class of dynamical r-matrices defined by Etingof and Varchenko 



In our paper we propose a very general framework for Hamiltonian reduction by pure 
gauge symmetries that allows one to derive a classical r-matrix. This framework is based on 
rather general assumptions: essentially that the first class constraints generate the adjoint 
action of some Lie algebra on a (Lie algebra valued) function on the unreduced space, which 
in turn reduces to the Lax equation, and that the Poisson brackets between these elements 
is already cast into r-matrix form. 

The organisation of the paper is as follows. In the second section we outline the general 
framework of the Hamiltonian reduction which allows us to derive r-matrices for integrable 
systems. In the third section we discuss the Dirac brackets and the r-matrices more gen- 
erally, relating each of these to generalized inverses. This sheds some new light on Dirac's 
brackets. The fourth section consists of our first example, where we apply our method to 
derive the classical r-matrix structure for the elliptic Calogero-Moser system with spin [21 [3 
using its Cech-like Hitchin description |lf)| [^j . The fifth section is devoted to the deriva- 
tion of the classical Feigin-Odesskii algebra [S| from some free field theory. This theory is 
a modification of the Higgs bundle related to the holomorphic GLAr-bundle of degree one 
over an elliptic curve. Instead of a Higgs field that takes values in the endomorphisms of 
the bundle we consider the field of automorphisms of the bundle. It is a Poisson manifold 
[0] and we show that the elliptic Belavin-Drinfeld r-matrix ^2] is naturally obtained in the 
context of such a reduction. In the concluding section we discuss a possible generalisation 
of our reduction technique to Hitchin systems on curves of higher genus. An Appendix is 
devoted to the special elliptic functions we use. 

Throughout the paper we use standard notations from quantum group theory and in- 
tegrable systems to express the Poisson brackets between entries of matrices. Thus, for 
example, the r-matrix equation cast as 



{L® L}= [r,L0l] - [r'^,lOL] 
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becomes 

{Li,L2} = [ri2,Li] - [r2i,L2]. (2) 

If the Lax matrix takes values in some Lie algebra g (or representation of this), the r- 
matrix takes values in g or its corresponding representation. We may expand quantities 
in terms of a basis of g as follows. Let denote such a basis with \t^^tu] = c^^ t\ defining 
the structure constants of g. Suppose 4'{t^) = X^, where (p yields the desired representation 
of the Lie algebra; we may take this to be a faithful representation. With L = L^X^ 
the left-hand side of (P) becomes 



{L^L} = ^{L^,L''}X^®X,, (3) 



while upon setting r = r^^X^ Xy and r = r^^^X^^ ® X^ the right-hand side yields 

[r, L ® 1] - [r^, 1 L] = r''^L^([X^, Xx] ® X, - X, ® [X^, Xx]) 

= {r^^4^L^-r^^4,L')X,(S>Xy. (4) 

Concretely, in the standard basis of gl]\f 
we have {X/^} = {^ij} and © takes the form 

m 

where 

f = ^ njkieij ®eki. (6) 

i,j,k,l 

Then r2i denotes the function with permuted indices r2i[z) = rkiij{z)eij ® e^i- 

At the outset let us record that under a gauge transformation r-matrices transform as 

+ iiAr^A^HAi, A2}, A2 '^2(^')A2] . 

Also, given a Lax operator, r-matrices are far from being uniquely defined. We will specify 
this ambiguity in due course. 



2 The Formal R-Matrix via Hamiltonian Reduction 

In this section we show that if a classical integrable system is obtained by a Hamiltonian 
reduction involving purely gauge symmetries, and so corresponding first class constraints, 
then the system admits a canonical r-matrix. The word "formal" in the title of the section 
stresses the fact that in many interesting cases the reduction is performed in the context 
of a field theory, and consequently for such systems the formula for the r-matrix must be 
properly defined. 

Our approach makes three quite general assumptions regarding the Hamiltonian reduc- 
tion that allows us to calculate the r-matrix. These are: 
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1. There exists an element which upon reduction becomes the Lax matrix, and this 
takes values in some Lie algebra g. 

2. The gauge transformations of <I> S g can be represented as 

$h^/i~i^>/i, (8) 

where h is assumed to take values in the Lie group G corresponding to the Lie algebra 
0- 

3. The Poisson brackets between the "entries" of the element $ in the unreduced space 
have been cast into the r-matrix form before reduction. 

The first two assumptions are quite natural if we wish to obtain an integrable system 
together with its Lax representation. The second assumption implies that the Poisson 
brackets between the first class constraints generating the gauge symmetries, and ^ can 
be written in terms of the commutators 

{Ta,<^>} = [ea,<^], (9) 

where {ea} is a basis of the Lie algebra of gauge transformations (which may be a subalgebra 
of g). The third assumption is more elaborate. It means that there exists a classical r- 
matrix G 0^0; defining the Poisson brackets in the total or unreduced space between 
the "entries" of $ 

{<Di,ci>2} = [r?2,^'i]-[r2°i,<I>2]. (10) 

Our first calculation below shows how the initial r-matrix produces, under reduction, 
the desired classical r-matrix of the integrable system. One may feel that we have simply 
transferred the problem of constructing an r-matrix for a reduced system to one of con- 
structing an r-matrix for an unreduced system. In practice constructing the latter r-matrix 
is often easier, but we will deal with the construction of in the next section. 

Together with some gauge fixing conditions = first class constraints form a 
system of the second class constraints {(Ja} = {Ta,x''} Ml- Using these constraints we 
may define the Dirac bracket, which allows us to calculate the reduced Poisson brackets in 
terms of the Poisson brackets on the unreduced phase space. Thus, in order to calculate 
the Poisson bracket between two observable /, k of the reduced phase space, we have first 
to calculate the Dirac bracket between arbitrary continuations F, K of f and k in the 
unreduced space: 

{F, K}db = {F^ K} - {F, ua} C^^{ap, K}. (11) 

Finally the reduced Poisson bracket between / and k is obtained by restricting expression 
Hll|l to the surface of the second class constraints {a a = 0}. Here ||C"^|| is the matrix 
inverse to ||Cq/3|| = {aa-,(yp} and { , } is the Poisson bracket on the unreduced space. 

Since under reduction the element $ € g becomes the Lax matrix, the Poisson brackets 
between the "entries" of the Lax matrix are just the on-shell Dirac brackets between the 
respective "coordinates" ^ on g. Our third assumption automatically casts the first term in 
the respective Dirac bracket (|11|) into the required form (|in|) . The perhaps surprising thing 
is that the remaining terms of the on-shell Dirac bracket Hll|) may also be cast into r-matrix 
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form as a consequence of the second assumption. To see this we observe that (on-shell) C 
takes the block form 

C-( ° {Ta,x'}\^( P \ pb-^T yH (12) 

using the fact that {Ta,Ti,} = as they are first class constraints. Then has the form 



{ P-' j = \P~' J 



(13) 



where 

Q''' = {p-'rc{x',x'KP-'t (14) 

Thus the only non- vanishing on-shell entries of the matrix ||C"'^|| are 

C>^"^i> = -C^ftX" = (p-i)b^ together with C^-^'^ = Q"^ . 
Using these the final term in {<I>i, ^2}db may be written in r-matrix form as follows: 
-{1>i, CT„} C^P{ap, $2} = -{1>i, xn(^-')am, ^2} + {^1. Ta}{P-^)t{x\ ^2} 

= -[{P-')lea ® {X\ <J>}, ^1] + [{X\ ^}{P-'t ® 65, f 2] 

Here we have made use of equation @. Upon combining all of the terms pO() for the 
on-shell Dirac bracket {$i,$2}db we obtain the desired classical r-matrix. 

Theorem 1 Under the three assumptions given we have 

r = (rO - {P-^Tbea ® {x', + W'^a ® {Tft, $}) (15) 

Z on shell 

where P is given in U^) and Q in Jj^l ). 

The terms modifying have the same general form as the infinitesimal form of We 
shall further explain this remark later in the paper. 

In section four we apply the method outlined here to derive the classical r-matrix for 
the elliptic Calogero-Moser system with spin. Before so doing however, let us return to the 
question of r^, and place our construction in a somewhat broader context. 



3 Dirac Brackets, Generalized Inverses and R-Matrices 

Our assumptions have a priori assumed the existence of an r-matrix r^ for the unreduced 
system. In many examples such unreduced r-matrices are easy to construct: when the 
unreduced matrices $ depend only on half of the phase space variables, such as the momenta, 
then for example r'' = is a possible solution. Such however is not always the case, and 
in this section we shall discuss the issue of the existence of r^. This will enable us both to 
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place our Dirac Bracket construction in a wider context and to elaborate on the remarks of 
|17j concerning this reduction procedure. For concreteness we shall phrase our discussion 
in the language of finite dimensional matrices (as we did in the previous section). 

The construction of r-matrices is essentially an algebraic operation. This algebraic 
nature was clarified in |1BJ where a necessary and sufficient condition was given for the 
existence of an r-matrix based upon the fundamental Poisson brackets of the Lax matrix 
(here The novel part of this investigation was the use of generalized inverses (and the 
construction thereof for generic elements of the adjoint representation of a Lie algebra g). 
We shall recall these notions in our present setting, relating them to Dirac brackets. In so 
doing we will arrive at some new results pertaining to Dirac brackets. 

We begin with generalized inverses. Let A be an arbitrary matrix. (In particular A 
need neither be square nor invertible.) A matrix Ga is said to be a generalized inverse of 
A provided 

AGaA = A. (16) 
Such a matrix always exists, yet need not be unique. We can further require 

GaAGa = Ga. (17) 

Again such a generalized inverse always exists, yet need not be unique. Observe that given a 
Ga satisfying (fTKI) and (|T7|) we have at hand projection operators Pi = GaA and P2 = A Ga 
which satisfy 

APi = P2A = A, PiGa = GaP2 = Ga- (18) 

One may specify a unique generalized inverse (which always exists), the Moore- Penrose 
inverse, by additionally requiring 

{A Ga)^ = AG A, (GaA)^ = GaA. (19) 

Let us remark that the adjoint ()^ here is defined with respect to a given inner product, 
the Moore-Penrose inverse satisfying a norm-minimising condition. Typically this is an 
hermitian inner product and the adjoint is the hermitian conjugate. We shall denote by 
A~^ the Moore-Penrose inverse of A. (Accounts of generalized inverses may be found in |191 
I2U1I211[^ .1 Geometrically the Moore-Penrose inverse may be constructed by orthogonally 
projecting onto the subspace on which A has maximal rank and inverting the resulting 
matrix. We also record an alternative characterisation |23j . Denote by A and G the matrices 

^ _ fO A\ n - ( ^ 
^ " Vo ' ^ " VGa 

Then 

1. (Uni) and jni) (^,(5, [i,(5]) form an s/2-triple. 

2. Further Ga = -4+ <^=^ [A, G] is Hermitian. 

Let us now relate Dirac brackets to generalized inverses. We may suppose our (say 
unreduced) phase space has canonical Poisson brackets 

{x\p,} = 5] (20) 
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so that 



{F,K} 



dF dF 
dx ' dp 



(dK_\ 

dx 



J 



J 



1 
-1 



(21) 



dK 
\ dp / 

(We shall give a coordinate independent description in due course.) In terms of the second 
class constraints we consider the matrices 



doa 
dpk ' 



dx^ ' 



Cap = {<7a, CT/j}. 



(22) 



In general we have a different number of constraints c from the number of coordinates, 
and the matrices a and b are not invertible. Given that the Moore-Penrose inverse is 
constructed by an orthogonal projection, and that Dirac's brackets give us a projection 
onto the constraint surface, the fohowing may not be surprising (though it should be more 
widely known). 



Theorem 2 With the definitions \2'J\) let A = (b, a) and set Ga 



(^) 

AGaA = A, GaAGa = Ga 
(ii) If Pi = GaA then Pi J = JP[. 
(Hi) 

'dF dF 



and 



AGA = ldc 



_^T^-i ) • Then 



(AGa)^ 



{F,K} 



DB 



dx ' dp 



(l-Pi)J 



dx 

dK 
\dp / 



(iv) Pi is self-adjoint with respect to any inner product of the form 

((u, v)) = u^Pf QPiv + u^(l - Pf )Q(1 - Pi)v 

where Q is an invertible symmetric matrix. For such an inner product Ga = A~^ is 
the Moore-Penrose inverse of A. 

The proof of these statements is straightforward. Observe that from the basic Poisson 
bracket 

C = ba^-ab^. (23) 



we have that 



ib,a) 



-b^c-^ 



ab^)C-^ = Ic 
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t ^ / a^C-' 



Thus P2 = AG A = P2 (foi' s-iiy choice of the adjoint). It is then clear that Ga = y 
satisfies ([TH]) and (|T7|) . Direct calculation shows that 



J, _n A-( ^^C-H a^C-^a \ 



whence Pi J = JPf , and {in) follows from (|1H) . We will have established the theorem upon 
showing Pi is self-adjoint. It is now that the issue of inner product confronts us. We shall 
discuss this more generally after proving the final assertion of the theorem . Now using the 
specified inner product and that Pi, P]- are projectors we have that 

((Piu, v)) = u^Pf • Pf QPiv = u^Pf QPiv = u^Pf gPi • Piv = ((u, Piv)). 

Then upon comparing with the definition of the adjoint, ((Piu,v)) = ((u, P|v)), we see 
that Pi = pI and so is self-adjoint. Therefore Ga = is the Moore-Penrose inverse for 
such an inner product. □ 

Let us now describe this theorem in a more geometric manner. Denote the unreduced 
phase space by (P, w). Given a function P on P we have a corresponding vector field Xp 
given via 

dF = ixpio = td{XF, ). 
The Poisson bracket of two functions is then 

{P, K} = u;{Xf, Xk) = dF{XK). 

Thus with uj = ^ dq' A dpj we have 

dF d dF d 

Xp 



dpj dxj dx^ dpj 

and the Poisson brackets H21|) . Our constraints describe a reduced phase space V = 
n'^^i{aa = 0} C P. Important for us is that V is not an arbitrary submanifold but a 
symplectic submanifold, and the inclusion i : V ^ P gives the symplectic form ajy = i*'^ on 
V. Because ujy is nondegenerate, TzV fl (T^l/)'^ = 0, where {TzVY = {u ^ TzP \ ijj{u,v) = 
G y}. Then 

r,p = TzV e {TzVy 

and we have a projection vr : TP — > TV . One finds that 

AXf) = Xf- {P, da} C"^ X,^ =Xf- X,^ dcj^{XF) = [1 - X,^ C7^" dcj^]{XF) 
and the Dirac brackets are 

{F,K}db = {F\v.K\v} = u{tt{Xf)MXk)) = dF{^{XK)) = dF{[l-X^^ C^^" 

(Observe that vr^ = vr follows from (i(TQ,(Xo-^) = Cap. Such a projection operator will exist 
even for a Poisson manifold P provided Cap is invertible: in this case we have T^P = 
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kervr Imvr and Imvr is again Poisson.) Comparison with the corresponding coordinate 
expression given by {in) of the theorem then shows that 

TT = 1 - Pi. 

The projection operator Pi is then Pi : TP (TzV)'^ , which is spanned by the vector fields 
{Xfj^}. Everything thus far has taken place within the realm of symplectic geometry: Pi is 
a symplectic projector and is self-adjoint with respect to the nondegenerate bilinear form 
u, 

uj{Pi{Xf),Xk) = uj{Pi{Xf),Pi{Xk)) = oj{Xf,Pi{Xk)) = u{Xf,pI{Xk)). 

(This is the content of {ii) of the theorem with u;(u, v) = u-^Jv.) To talk of orthogonal 
projection we need an inner product. While the kinetic energy of a natural Hamiltonian 
system can provide this an inner product is an additional ingredient. However, as {iv) of the 
theorem shows, there are a large class of inner products for which tt becomes an orthogonal 
projection and we have 

T,P = T,V®{T,V)^. (24) 

In practice one often further restricts attention to compatible inner products for which we 
have 

((u, v)) = a;( Ju, v), lj(u, v) = a;( Ju, Jv), = — Id 

and respecting the vector space decomposition (^l)) . 

We remark that generalized inverses have been discussed in various connections with 
singular Lagrangian systems. Broadbridge and Petersen f24_ refine an observation of Duffin 
|25j that a generalized inverse may be used to go from a singular Lagrangian system, such 
as those arising when constraints are implemented by Lagrange multipliers, and a corre- 
sponding Hamiltonian system. The reduced Hamiltonians appearing correspond to those 
of de Leeuw at al Q^. (These works place various restrictions on the nature of the con- 
straints.) Because the Moore-Penrose inverse may be calculated very efficiently their use 
in solving constrained dynamical systems is important. Kalaba and Udwadia formulate a 
large class of constrained Lagrangian systems as a quadratic programming problem, and 
use the Moore-Penrose inverse to solve these |27) . 

Finally, just as the Dirac bracket can be expressed in terms of generalized inverses, so 
too can the solution of the r-matrix equations (|2j or (|lflj) . From Q the r-matrix equation 
takes the explicit form 

b = a^r — r^a 

where we have set a^^ = c^^xL^ = —ad{LY^ and h^^ = {L'^,L'^}. The solutions of this 
equation have been studied 28 . 

Theorem 3 Let g be a generalized inverse of a = —ad{L) satisfying U6}) and jlTj ). Then, 
with Pi = ga and P2 = ag, the r-matrix equation ^\) has solutions if and only if 

(1-Pf)6(l-Pi) = 0, (25) 

in which case the general solution is 

r = ^g^bPi + /6(1 - Pi) + (1 - Pj)y + {P^ZP2)a (26) 

where Y is arbitrary and Z is only constrained by the requirement that P2ZP2 be symmetric. 
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The theorem then constructs the r-matrix in terms of a generahzed inverse to adL and 
describes the ambiguity of the solution. Although the general solution appears to depend 
on the generalized inverse, the work cited shows that changing the generalized inverse only 
changes the solution by such terms. Further, the generalized inverse of a generic adL may be 
constructed. Thus the existence of an r-matrix has been reduced to the single consistency 
equation ()25|) and the construction of a generalized inverse to ad{L). 

Having now addressed the issue of , and indeed the general construction of r-matrices, 
we will return to the particular construction given by the Dirac reduction procedure of the 
previous section. We shall illustrate this in turn by two examples. 

4 Example 1: Derivation of the R-Matrix for the ElUptic 
Calogero-Moser System with Spin 

In this section we show that a classical r-matrix for the elliptic Calogero-Moser system with 
spin can easily be derived with the help of Dirac brackets. One surprise (for us) is that in 
calculating the r-matrix we do not make a use of the identities for elliptic functions that 
typically under ly these integrable systems. 

In deriving the r-matrix for the elliptic Calogero-Moser system with spin we will divide 
the Hamiltonian reduction procedure into two stages. In the first stage we omit gauge 
transformations related to the parabolic subgroup. This yields the formula for the classical 
r-matrix of the spin Calogero-Moser system with the maximal spin sector. 

After completing the first stage of the reduction we are left with residual action of 
the parabolic subgroup on a finite dimensional phase space. Here it is possible to derive 
the desired classical r-matrix using any concrete fixing of the residual gauge symmetry. 
An alternative method also exists, making use of a gauge invariant extension of the Lax 
matrix. We will adopt the latter approach to show that upon further reduction the r-matrix 
obtained in the first stage becomes a classical r-matrix for the elliptic spin Calogero-Moser 
system. 

The Hamiltonian reduction that leads to the elliptic Calogero-Moser system with spin 
was originally presented in the paper [S], where the phase space of the system was realised 
as the cotangent bundle to the moduli space of topologically trivial holomorphic bundles 
(the Higgs bundles) over the torus with a marked point. We are going to employ this 
construction in deriving the respective classical r-matrix. 

For our purposes we mostly follow the paper W. and use a Cech like description of the 
moduli space. We will realise the torus S,- as a quotient of C/gZ with q = e^'^*'^, Imr > 0. 
We choose as the fundamental domain the annulus Ann-r = {\q^^'^\ < \z\ < \q~^^'^\}- 

We define a holomorphic vector bundle En over S,- by the matrix transition function 
g{z) G GLn{C) that is holomorphic in a neighbourhood of the contour 7 = {\z\ = \q^^'^\}- 
In this way g{z) represents a Cech cocycle. The set of these fields we denote as C = {g}. 

Let G be the gauge group of holomorphic maps f{z) of the annulus Anrir to GL]\j[C) 
such that 

f{z) , = /, (27) 
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where I is the identity matrix. It acts on C as 

g{z)^ f{z)g{z)r\q-^z). (28) 

Observe that the Lie algebra Lie(t/) consists of matrix-valued holomorphic functions £{z) : 
Anrir I— > gl]\j{C) vanishing at the point z = 1 



= 0. 

2 = 1 



e{z) 

The quotient space under the gauge action (|28() 

M = C/g (29) 

is the moduli space of the holomorphic bundles over E,- with the marked point z = 1. Note 
that upon reduction with respect to adjoint action of a parabolic subgroup V C SL]\f{C) 
the space H29() becomes the moduli space of the holomorphic bundles over with a quasi- 
parabolic structure at z = 1. 

The cotangent bundle T*C is called the Higgs bundle. The dual field (the Higgs field) 
is a one-form ^ 

V = viz) — 

taking values in End (-Eat) = Lie algebra gl]\f. The field rj{z) is holomorphic in a neigh- 
bourhood of 7i. 

The bundle T*C is the non-reduced phase space with the symplectic form 

= ^ jf tr5(g-i (z)r?) A 5g{z) . (30) 

We lift the gauge group action (PH|) on T*C as 

7^iz)^f{zUz)f-\z). (31) 

The symplectic form ()3U() is invariant with respect to this gauge action. In this way it 
produces the moment map 

fi : T*C^ Lie*{g), (32) 

which is defined as 

f^[Tl,9]{^) = f i^{9"\z)V9{z)e{q'^z) -rie{z)), e ^ Lie{g). (33) 
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The zero level of the moment map ()33p is the surface in T*C where rj[z) can be extended 
meromorphically from the boundary (71, 72) inside the annulus with at most single pole of 
the first order at the point z = 1 and such that for \z\ = \q^/'^\ 

rjiq-^z)=g-\zHz)g{z). (34) 

While a generic topologically trivial bundle over a torus can be decomposed into a direct 
sum of line bundles |29: , a generic holomorphic bundle with a marked point can differ from 
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a direct sum of line bundles by the conjugation of a constant matrix a S SLi\i{C). In 
other words, using the gauge group action ((2 
constant matrix of the following form 



/ g-27rmi 





a 





-2TTiU2 



one can transform a generic field g{z) to a 



v 



In particular, this means that 










-27viuj\ 



a 



I 



dzgaiz) 







(35) 



(36) 



is a proper gauge condition for our system. 

It is easy to see that the entries of the matrix a are defined by h up to the scaling 
transformations 



aij ^ OiijXj, Xj / 0. (37) 

while Ui and = Uj + ni + n2T , ni, n2 € Z gives the same holomorphic bundle over S,-. 
Thus, an open dense set of the moduli space M. (|29|) can be parameterised by the points of 
the set^ 

(J X CP^-i)^, (38) 

where J is the Jacobian of the curve S,-. In the general context of holomorphic vector 
bundles over Riemann surfaces this is known as the Tyurin parameterisation 30^. This 
parameterisation has found recent application in a general description of Hitchin systems 
and their r-matrices jSJ; ISS, |33] . 

Expanding the functions g{z) and r]{z) in a neighbourhood of the contour 7, 

rii.i^) = E ^^'^-^^ ' = E gtjz'^ (39) 

aeZ aeZ 

we readily find the respective Poisson brackets to be 

Hj,Vm} = -Vti^'h, + Vtt%i, {9?j,Vm} = gtt'^u, {9?j,9m} = 0. (40) 
In these terms the expansion coefficients of the constraint f^['r],g] = (|33)) are of the form 

^ij = E i9-'tk'vli9f-V - E (a-'tk^'viigr/ - Vtj + ^° , a / 0, (41) 

b,c,k,l b,c,k,l 

and the gauge fixing conditions H36|) can be rewritten as follows 

gfj = 0, a / 0. (42) 
Our second class constraints are then 

Mg. = ( g / 0), ff«=0(a^O). (43) 

^In fact we slightly abuse the parameterisation omitting the quotient with respect to the symmetric 
group of permutations Sn, which plays an analogous role as the Weyl group for the rational Calogero-Moser 
system. 



12 



The field g{z) being restricted to the constraint surface ()43() becomes the constant matrix 
g{z) = h (|35j) while the field r] yields the Lax matrix of the elliptic Calogero-Moser system 
with maximal spin sector, 

r]{z) = l{z) = al{z) a-^ (44) 



with 

N 

l..f-y\ — _ 

27rz' 



k=l 

Here Ui, Oij, vj and (3ij are coordinates on the reduced phase space T*M. ~ T*{J x 



QpAf i^AT^ Because aij (for each fixed j = 1, . . . , A'^) are homogeneous coordinates on 



N N 



the respective projective spaces CP , the symplectic form on T*JV[ is obtained from the 
form 



r^i = dvi A dui + dPij A daji (45) 

i=l i,j=l 

by symplectic reduction on the first class constraint surface (arising from the generators of 

N 

X]Afc«fci = 0- (46) 

k=l 

The function (f){z, u) in (|44() can be represented in the following form 

<P{z,u) = Y, aj.iu_i ^ |g|<|z|<l, u^a + hT, a,heZ. (47) 

aeZ ^ 

We note that the Lax matrix of the Hitchin system with a quasi-parabolic structure 
at the marked point z = 1 is obtained from 1)441) by reduction with respect to the adjoint 
action of the corresponding parabolic subgroup V C SLn{C) . 

In order to construct the r-matrix using our theorem we must determine the Dirac 
bracket. First we calculate the matrix of the Poisson brackets between the constraints ()43|) 
and then determine the inverse matrix. The non-vanishing on-shell brackets are found to 
be 

{4 , } = 6a+b,o {q'^hiidkj - hkj6ii ) . (48) 

To present the inverse matrix we adopt the following convention for the entries of the inverse 
matrix: (^(j^ay^ stands for the entry with Mf- being the first index and g^.^ being the 
second. Then the non-vanishing entries of the inverse matrix can be written in the following 
form 

((^Mg\ij kl _ _,^gM\kl ij _ ST^ ^a+b,0 Ctlmj'^ ^)miQjra(Q ^)nk 
V'-' )a b ~ )b a ~ a g-27rm„ _ g-27rjMm ' ^ ) 

m,n=l 

In deriving the r-matrix via the Dirac bracket we will exploit the following simple 
property 

{M-,ri'} = [e,,,r+']-[e,,,rj% (50) 
which is just the counterpart of the general relation ©. 
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Now the expression for the Dirac bracket between the expansion coefficients of the 
field r]{z) consists of two parts. The first part is just the initial Poisson bracket in the 
extended phase space. This may be rewritten as 

{^1,4} = \ En ® ^^-^ ^1^'] - \ E[«^^- ® ^3^^^t\ (51) 

The second part reflects the reduction. With {f?i,??2}£'S = {^^\^4l\ + the on-shell 
expression for the second part is 

m,n=l 
m,n=l 

where Umn = Um - Un- 

Finally, to calculate the Dirac bracket between the fields r]{z) and r]{z') we must to 
perform the sum 

{rj^{z),V2{z')}DB = E W,r?^}DB^"(/)^ (53) 

a.feGZ 

Such sums require some care because separately certain series in the expression for the 
Dirac bracket (|51I52|1 do not converge. These must be accurately combined in order to get 
a finite answer. 

For example, if we choose \q\ < \z'\ < \z\ < 1, hence \q\ < \z' j z\ < 1 and 1 < \z/z'\ < 
\q~^\, then the second sum in (|52|) gives a divergent series in (|53j) . However if we add to 
the above terms the expression ()51() we obtain a convergent series and the desired r-matrix 
for the Lax matrix (|44j) takes the following form 

r{z, z') = axa2{{E{z' I z) - E{z')) E en eii+ 

i 

+ '^{cl){z' /z,Uji) - (j){z',Uji))eij eji]a^^a2^ , (54) 

where 

ai = a(8)/, a2 = I ® a 
and the function E{z) can be represented in the form of the following series 

^(^)= E ^-1' H<I^I<1- (55) 

Remark 1. The classical r-matrices for the spin Calogero-Moser systems were originally 
found in the paper j^. In contrast to |34| where the final formulae, both for the Lax matrix 
and for the r-matrix, were defined on the elliptic curve only after an additional auxiliary 
reduction, our formulae 1)44(1 and (|54|) are defined on the elliptic curve from the outset. 
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Remark 2. The elliptic Calogero-Moser system with maximal spin sector can be also 
described in terms of a Lax matrix and r-matrix which are functions, rather than sections 
of some bundles on the elliptic curve. Namely, the Lax matrix ()44|) and the r-matrix H54() 
are gauge equivalent to the following matrix-valued functions on T,r 

N 

= ^ aikBkiiw){a~^)ij, Bii = 



k,l=l 

N 



(56) 

ft^^{w)9n{w - Ui)9n{ui - Uj)9n{uj) 



k=l 



N 

r{w, w') = ^ {E{w - w') + E{w'))eij ® eji~ 

N 



(57) 

^ aii{a~'^)ik{E{w - ui) + E{ui))eij ® ejk, 

i,j,k,l=l 



where 



w = In z , w' = In z' . 

2m 2m 



and 

9ii{w) = ^ exp {mT{m + l/2f + 2m{m + l/2){w + 1/2)). 

mGZ 

The Lax matrix for the spin Calogero-Moser system in the form H56|) was originally obtained 
in the paper by Krichever |S2], while the formula (|57|) for the corresponding r-matrix has 
been given recently in j33j . 

Remark 3. In obtaining the result ()54|) we have only used the definition of the Dirac 
bracket and not the elliptic function identities which are unavoidable for checking that (|54|) 
is indeed the r-matrix for the Lax matrix (|44() . Of course the functions our procedure yields 
as series are, as shown in the Appendix, expressible as elliptic functions. 

The resulting r-matrix (|54|) is universal in a sense that it allows us to derive a classical 
r-matrix for the respective Hitchin system on the elliptic curve with an arbitrary quasi- 
parabolic structure at the marked point. We will conclude this section by elaborating upon 
this point. 

Let V he a parabolic subgroup of SLn{C)- Then the Hitchin system on the elliptic curve 
with the corresponding quasi-parabolic structure at the marked point is obtained from the 
maximal spin system described above upon the further Hamiltonian reduction with respect 
to the following action of the parabolic subgroup V, 

N N 
o^ij ^'^QikCtkj , ctij ^'^Pik{q~^)kj , Ui^Ui, Vi^Vi. (58) 

k=l k=l 
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Here q = {qij) € V. The momentum map of the action 1)58(1 is found to be 

: T*M ^ Ue*{V) , fi^[a,f3]{X) = tv{apX) , X G Lie(P) C g/jv(C). (59) 

We note that the Lax matrix H44() transforms under the action of g E "P as 

l{z)^ll{z)=ql{z)q-\ (60) 

and this agrees with our general assumption Q. 

This stage of Hamiltonian reduction is finite dimensional and, whatever parabolic sub- 
group V and gauge fixing conditions are chosen, our formula (|15() yields the desired r-matrix. 
To proceed further with this final Hamiltonian reduction we must now choose a concrete 
parabolic subgroup V and impose the constraints fi^{a,f3) = together with some gauge 
fixing conditions. One could implement this simply using the Dirac reduction procedure 
once more to obtain finally an r-matrix for the elliptic spin Calogero-Moser system with 
the chosen quasi-parabolic structure. An alternative method however exists, and we will 
use this. This alternate route proceeds by constructing a gauge invariant extension of the 
Lax matrix. What is this gauge invariant extension? We have said that our Lax matrix 
transforms as ()6fl|) . A gauge invariant extension of this is constructed by conjugating by a 
compensating gauge transformation that "undoes" this gauge transformation. We imposed 
some gauge-fixing constraints x"^ = so as to choose (locally) one representative from each 
gauge orbit. Thus points in our total space can be described (locally) by coordinates {p,q), 
where p is a point on the constraint surface = Q describes the gauge transformation 
along the orbit. The gauge invariant extension of the Lax / is 

l^{z)=q-H{z)q 

where if l\p\ transforms to q~^l q at the point {p, q) then q{p, q) = q~^ so that l^{p, q) = l{p)- 
Therefore the twisted Lax matrix is a gauge invariant extension. Now we have described 
in Q how the r-matrix transforms under a gauge transformation, the desired r-matrix for 
the considered Hitchin system with the quasi-parabolic structure V are found to be 

r^{z,z') = {qi^q2^r{z,z')qiq2 

+ qi^q2^{H^')^Qi}Q'2 + l[qi\2^Ui^Q'i}^Q2^Hz')q2]) ^- (61) 

/ on shell 

Observing that q{p) = Id, this expression reduces to p5j) . and so explains the origin of the 
remark we made in deriving our theorem. Thus by calculating ()61|) via the gauge invariant 
extension we will obtain the same r-matrix as determined by Dirac reduction. Reduction 
to the spin zero Calogero-Moser models yields the r-matrices of [HSl IHBl 07] . 



5 Example 2: The Feigin-Odesskii Bracket Via Poisson Re- 
duction 

In a recent paper j^-iSj the classical Sklyanin algebra has been derived by Hamiltonian re- 
duction to the moduli space of complex structures on a principal GL2-bundle of degree one 
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over an elliptic curve. In this section we generalise the construction of [3H to the case of 
GLn in order to obtain the classical Feigin-Odesskii algebra 

Let be an elliptic curve and let En be a principal GLjy-bundle over S,- defined by 
the following gluing rules for a section s: 



s{z + 1) = hsiz), s{z + t) = A{z)s{z). 



(62) 



Here 



h 



(I 

e 



V 







A{z) = I2 exp 



-7V-1 



2'iriz\ 





1 






1 



\ 







V 1 








1 

/ 



and e = exp (^) . 

The Cech cocycle condition determining this bundle follows from the commutation re- 
lation 

hh = e-'hh. 

It easy to check that the holomorphic section of the corresponding determinant bundle 
det(£'Ar) is just a ^-function with a single simple zero on E,-. We therefore conclude that 
deg(£'Ar) = 1. Let us also note that 



la 



jai 70,2 
^1 ^2 



(01,02) e Ziv X z 



N 



(63) 



forms a basis oi gljy. 

In our construction of the Feigin-Odesskii bracket we use deformations of the complex 
structure on that preserve the determinant bundle det(i?7v) of E^. In contrast with 
previous Section we will now use the Dolbeault picture to describe the complex structures 
on En- The operator 



d^ = kd + A : ^P''?(S^,Sjv) ^^^•^+HSr,^jv) 



trA{z,z) = 



(64) 



acts on the sections of En and defines a complex structure on En- Here A = A(z,z)dz is 
a (0, l)-connection of En- 



A{z + 1) = hA{z)l^ 



A{z + t) = hA{z)L 



-1 



(65) 



The constant k appearing in (|64j) may be further identified with a central charge. Note 
that the traceless condition trA{z, z) = guarantees that the determinant of the deformed 
holomorphic bundle En, defined hy kd + A coincides with that of En- 

Two complex structures A and A^ are called equivalent if they are related by the fol- 
lowing transformation 

A^Af = f-'Af + f-'kdf (66) 
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where / = f{z,z) is a smooth SLTv-valued function on St- satisfying the following quasi- 
periodicity conditions 

f{z + 1) = Iif{z)I^\ f{z + t)= l2f{z)I,\ (67) 

Now we define the total Poisson space P that will be reduced as a principal affinization 
over a cotangent bundle to the space of smooth sections of Aut{E]\}). A point on the 
respective base is identified with a GLjy-valued field g = g{z,z), satisfying the following 
quasi-periodicity conditions 

g{z + 1) = Iig{z)I^\ g{z + r) = hgiz)!^^; (68) 

and a point on the fibre of P is determined by a (1, 0)-connection A ()65|). 
The Poisson brackets on P are defined in the following way: 

{Aij{z,z),Aki{w,w)} = {Au{z,z)6kj - Akj{z,z)6ii) 6{z - w) 5{z - id) 

+ k{6ii5kj - -^Sij6ki)d6{z - w) 6{z - w), 

{gij{z,z),AM{w,w)} = {gii{z)6kj - ^gij{z)5ki) d{z - w)5{z - w), (69) 
{gij{z,z),gki{w,w)} = 0, 

where A; is a central charge. The space P with this Poisson structure is a particular case of 
a general construction proposed by Polishchuk [S]. If we supplement the transformations of 
A ()66() with the simultaneous transformation of the field g 

A^Af = f-^Af + f-^kdf g^gf = r'gf (70) 

then the equations (|7Ujl define a Poisson action on P. 

Let us consider the Poisson reduction with respect to the action ((TOJ. First we note that 
a generic field A is in fact a pure gauge^ 

A = r^kdf (71) 

and we can choose ^ = as an appropriate gauge condition. Further, there are no residual 
gauge symmetries. To see this, consider a holomorphic function / : C i-^ SL]\f satisfying 
the quasi-periodicity conditions l|E7|) . By expanding in the basis l\&6\i we find that 

faiz + 1) = £-"2 faiz), fa{z + t) = s'^^ /„(z). 

Such a function must be a constant f{z) = Id. Thus the desired Poisson quotient is 
parameterised by a smooth section L = L{z, z) of Aut{Ei\i) 

L{z,z) = f{z,z)g{z,z)r\z,z). (72) 

In order to calculate Poisson brackets on the reduced space we have to construct a gauge 
invariant functional L = L[g, A] which coincides with L{z, z) on the surface of gauge fixing 

^In other words this means that there are no moduli of (semi)stable complex structures on En- 
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^ = 0. For this we solve the equation (|71|) for f{z, z) and substitute the solution into (|72() . 
Thus we get the desired gauge invariant extension 

Lb, A]{z, z) = f[A]iz, z)g{z, z)r\A]{z, z). (73) 

The calculation of the Poisson bracket between the entries of the matrix (|72|) is reduced 
to the calculation of an on-shell expression using the Poisson bracket ()69() between the 
entries of ((75|) . Note, that in doing these calculations we do not need to know the explicit 
dependence of on A. The only terms that enter the on-shell expression are 



rijkiiz,z;w,w) = k 



^fij{z,z) 



6Aik{w,w) 



-k 



6{f-%{z,z) 



6Aik(w, w) 



(74) 

A=0 



Due to equation (|7T1) the function rijki{z, z; w, w) turns out to be a Green's function for 
the (9-operator 

dzrijki{z,z;w,w) = {6ii6kj - ^dij6ki)d{z - w)6{z - w). (75) 

Note that the function 1)741) can be found in the form r(z, w) = r{z — w) and "unitarity" 
in the sense that 

r{z) = -r2i{-z). (76) 

In virtue of equations (|65j) . ()67)) . we have that r{z) also possesses the following quasi- 
periodicity conditions 

r{z + l) = {h(g)Id)r{z){I{^ 0ld), r{z + T) = {l20ld)r{z){I^^ (^Id). (77) 

Thus rijki{z) defines a meromorphic function r{z) : C i— > sIn (8) sIn with at most simple 
poles only at the points z = n + mr n,m £ Z. The residue t at the point z = is the 
Killing form of the Lie algebra si m 

t = '^{eij iS> eji - ^eii0 ejj). (78) 



Following the paper [T^ there is a unique meromorphic function r(z) satisfying the above 
conditions and, in particular, it is a solution of a classical Yang-Baxter equation 

[ruizi - Z2),r23{z2 - zs)] + [ruizi - Z2),ris{zi - Z3)] + [rn{zi - Z3),r23{z2 - Z3)] = 0. (79) 

Returning then to the calculation of the on-shell expression for the Poisson bracket 

{Li[g, A]{z, z), Lab, M^, w)}\a=o 

there are two types of contribution. The first type of contribution originates from the 
Poisson bracket between g{z) and /[^](t(7), and the respective expression is 

Bi = 2k~^Li{z, z)L2{w, 'w)r{z — w) — k~^Li{z, z)r{z — w)L2{w, iv) 

-k-'^L2{w,w)r{z -w)Li{z,z). (80) 



19 



The second set of terms arise from the Poisson bracket between f[A](z) and f[A]{w). These 
have the fohowing form 

B2 = —k~^r{z — w)Li{z, z)L2{w, w) + k~^Li{z, z)r{z — w)L2{w, w) 

+ k~^L2{w, w)r{z — w)Li{z, z) — k~^Li{z, z)L2{w, 'w)r{z — w). (81) 

Combining the expressions (|8U() and (|81|) we obtain the classical quadratic r-matrix 

{Li{z, z),L2{w, w)} = z) (8) L{w, €i),r{z - w)] (82) 



The Jacobi identity for the Poisson bracket (|82j) follows from the classical Yang-Baxter 
equation ((7^1) . 

In order to get an explicit expression for the r-matrix (|74|) we use the basis (|63|) and 
certain automorphic functions on C satisfying the following quasi-periodicity properties 

ifa{z + l)=e-''^^a{z), ^a{z + T)=e''^^a{z). (83) 

It is easy to see that up to an overall factor for each a ^ there exists a single function 
^Pa satisfying the above periodicity conditions and having at most simple poles only at 
the points n -|- mr, n, m € Z. For a = (0,0) we simply set '^{Qfi)[z) = 1. Then the elliptic 
r-matrix takes the following form 

^(^) = ^ Yl e^^^^^a{z)IaXl-a. (84) 

aeZjvxZjv 117^0 

Let us return to the Poisson bracket (|82|) . Consider the following set of functionals 

Cij= [ d^z i;{z,z)dLij{z,z) i,j = l,...N, (85) 

for all smooth functions ip vanishing at the points n + mr, n,m (z Z. We will prove that 
this is a set of weak Casimir functions^ for the Poisson bracket H82() . 

To this end let C be the matrix with entries H85() . Straightforward calculation then 
shows that 

{C l,L2{w,w)} = / (fz 'ip{z, z)d{Li{z, z), L2{w,iB)} 

= I Sziljiz,z)d\\Liz,z)(^L(w,w),r(z-w)]. (86) 

Up to terms vanishing on the surface Cij = (for all i, j) the expression (|86|) takes the 
form 

{C (g)l,L2(w,id)} = I d^z i;(z,z)hL(z,z)(^L(z,z),t6(z -w)6(z -w)] (87) 



^Functions Ca are said to form a set of weak Casimir functions if they generate a proper ideal in the Lie 
algebra of smooth functions with respect to the Poisson bracket. 
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where t is the residue of the r-matrix r{z) at the point z = 0. Since t is adjoint invariant 
we see {C 1, L2{w, w)} = 0. A consequence of this result is that the zero level surface for 
the functionals (jH^ 



Cij = 

is Poisson. That is, one can naturally define a Poisson bracket on the surface (|88() such that 
the embedding of the surface into the initial space is a Poisson map. 

Now the surface ()88() is finite dimensional. In particular, the matrix L{z, z) satisfying 
the equations (|88|l can be represented as 

L{z)= Yl SaMz)Ia, (89) 

where Sa are c- numbers parameterising the surface (|88|) . Using some obvious properties of 
the automorphic functions ifa one can easily show that the relation H82|) defines a quadratic 
Poisson bracket between the coordinates Sa ■ If we substitute the solution (|89|) in equation 
(|82jl this is just the desired Feigin-Odesskii bracket, which is defined as a classical limit of 
the Feigin-Odesskii algebra 0- 

Let us also note that all of the Casimir functions of the Poisson brackets (|82|) , (|89j) can 
be constructed with the help of our proposed reduction procedure. To see this we note the 
functional det g{z) yields a continuous set of Casimir functions for the Poisson brackets (|69|) 
of the unreduced space. Because the functional det g{z) is also gauge invariant it defines 
a continuous set of Casimir functions det L(2:) for the Feigin-Odesskii brackets (|82|). (|89j) . 
The desired algebraically independent Casimir functions for 1)821) . (|89|) may then be defined 
as the following coefficients of the Laurent expansion for det L{z) around the marked point 

Co = y"^^^^, C-2 = j zdetL{z), ... C-N = j z^-^ det L{z), (90) 

where F is a small contour around the point z = 

To conclude this section we would like to mention a relation between the matrix H89() 
and the Poisson brackets (|82|) . and the higher dimensional elliptic top j^lj. The equations 
of motion for this integrable system are usually given in the form 

|L(.) = i{L(.).C-.}|,, (91) 
where {•, ■} stand for the following linear Poisson brackets 

{Li{z),L2{w)}\^ = i([Li(z),r(z - w)] - [L2{w),r2i{w - z)]), (92) 

and C-2 is a quadratic Casimir function (|9U() of the Poisson brackets 1)82(1 . ()89() . The same 
equations of motion ()91j) can however be written in an alternative manner with the aid of 
the quadratic Poisson brackets H82() : 

f^L{z) = {L{z),H} (93) 



It is easy to see that the residue of the function det L{z) at the point z = is vanishing. 
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The Hamiltonian now has a very simple form 

H = —Sofi- 

Note that the natural generalisations of the elliptic tops [3^] associated to elliptic curves 
with multiple marked points can be obtained analogously in the framework of the our 
construction. For this we must require the functions 'ip{z,z), defining the weak Casimir 
functions (|85|) vanish at each marked point. 

6 Concluding Remarks 

In this paper we have shown how the technique of Hamiltonian reduction enables us to 
calculate rather than guess some important ingredients in the theory of integrable systems. 
Although the technique is not always rigorous in the field theory context where conditionally 
convergent sums frequently arise, this approach does give us some important geometric 
information about the objects in question. 

A stumbling block for the generalisation of our results to curves of higher genus is the 
lack of a convenient basis both for meromorphic functions and for meromorphic differentials 
on a general Riemann surface. Thus in the paper where a classical r-matrix is presented 
for Hitchin systems on an arbitrary Riemann surface of genus (7 > 2, it is mentioned that a 
formal expression for the r-matrix can be given in terms of a series in the Krichever-Novikov 
type basis 00]; While we have shown that sums of the type (jHSj) can be performed 

for the elliptic curve, for an arbitrary curve of genus g > 2 the series yields only a small 
amount of information about the geometric nature of the r-matrix in question. 
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Appendix: Representations of the Special Functions (j){z, u, q) 

and E{z, q). 

The meromorphic function u) is defined by the following properties. First, it is auto- 
morphic with respect to the transformation z 1— > qz, namely 

(/)(g-iz,n) = e2™0(z,u). (94) 

Second, on the unit circle |2:| = 1 the function has one simple pole at the point z = 1 and 
the respective residue equals 1. 
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It is easy to check that the function with the above properties can be represented as 

" 2^ + ^^""^ " + f ^^^^""^ " " ^^""^^ + •••)• ^^^^ 

Here 

(9^^(t^, t) = Y^ exp(7rir(m + 1/2)^ + 2m{m + l/2){w + 1/2)) (97) 

meZ 

and z = e^'^*"'. The function is the unique odd theta function, u, and p are the Weier- 
strass functions respectively of the same names, and rji = C(^)- Due to the automorphic 
properties of the 6'-function H97|) the right hand side of the equation is in fact a function 
of z and q = e^'^^'^ . 

On the other hand, using the above properties one can find series representations of the 
function (|95|) and thus prove that the properties uniquely define the desired meromorphic 
function. Due to the poles of the function 1)95^ on each of the circles \z\ = \q\"^, m G Z the 
series representations depend on the choice of the annulus. In deriving the r-matrix we use 
two annuli, 1 < [zj < \q^'^\ and \q\ < \z\ < 1. The series representations of ()95() in these 
annuli are of the form 

<Piz,u) = a^-i^iu _ 1 ' l<kl<k"^l> u/a + 6T, a,6GZ, (98) 

and 

(t){z, u) = Y ^ ^^.^ _ ^ , |g|<|z|<l, u^a + 6r, a,6G Z. (99) 

aeZ ^ ^ 

The meromorphic function Ei^z^ is defined by analogous properties. It is automorphic 
with respect to the transformation z i— s- q~^z 

E{q~^z) = E{z) + 1 (100) 

On the unit circle |z| = 1 the function has one simple pole at the point z = 1 and the 
respective residue equals 1. Finally, the function E{z) has a vanishing moment in the 
annulus 1 < \z\ < \q~^\, namely 

f]z 

-—E{z) = 0. (101) 

|z| = |g|-l/2 2-KlZ 

Using the ^-function (|97jl we may represent E{z) as 

E{z^ q) = ^^4^ -\ = ^ (C(^) - 2m - i. (102) 

The series representations for E(z) in the annuli 1 < [z| < \q^^\ and \q\ < |z| < 1 are of 
the form 

^(^)= E i<N<r'i, (103) 

aeZa^O ^ 
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and 
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